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I. INTRODUCTION 



There is an assertion that classical electromagnetic fields can be used to control the 
quantum behavior-i^. It has the justification from recent experimental works showing that 
classical fields can be used to excite the population of atorn^ or to extremely accelerate the 
neutral atoms^. In this paper, we show that the generation of quantum gates is necessary due 
to the presence of classical electromagnetic fields on a moving n-relativistic spin quantum 
bits. The proposal leads to the possibility to open the definition of quantum transistor since 
the generation of quantum gate is related to the certain bandgap opening. 

Bandgap engineering is playing crucial role in classical information theory-. To generate 
the classical bits in a transistor, the potential should be tuned by step functions determining 
the conductivity of the gate in transistor: if the potential exceeds threshold of bandgap, the 
transistor in conductive state meaning generates "ON" state, while if it is less than the 
bandgap threshold, the transistor is in a non-conducting state representing "OFF" state. 
Conventionally, these two states are implemented by constant functions. Here, we present 
that the notion of Dirac four-potential can be exploited to emerge the concept of band gap 
in a quantum transistor: a device used to generate quantum gates due to the presence of 
electromagnetic fields on the relativistic spin qubit in which the perturbation is represented 
by Lorentz force on the Dirac Hamiltonian. Contrarily with the classical transistor, the 
proposed quantum transistor has dependencies of the space-time structures in the Dirac 
potential. 

Following Moore's Law, the present efforts attempt to minimize the size and increase 
the capabilities of transistor. However, physical complications, for instance quantum effects 
such as the Kondo effect^ and quantum interference^, are increasing due to the miniatur- 
ization. Several theoretical and experimental frontier efforts have been proposed, in ad- 
dition to single-electron transistor-, lateral quantum dot^, single-photon transistor—, and 
single-molecule optical transistor—, nevertheless the science and technology for a transistor 
performing quantum computation is still unclear— . 

The strength of quantum computation for solving problem of computation have been 
known in the recent years^. Recently, there are several attempts to establish an appropriate 
theory to encode qubit and to implement quantum computation in realistic physical systems 
along with measurement-based quantum computingi^ and adiabatic quantum computingi^. 
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However, theory of how to implement it into the physical system involving the concept 
of quantum transistor considering the relativistic effect is still poorly comprehended. To 
date, we only find that only Ref.— is related to quantum transistor in which can potentially 
implement quantum computation operation for one qubit. 

Recent works on quantum computation are constructed based upon Schrodinger Hamil- 
tonian at which the potential is scalar. Thus, the perturbations of the potential due to 
the presence of external forces, such as Zeeman and Stark effects, are expressed in the 
scalar form. Lately, following the successful application of two-dimensional Dirac equa- 
tion describing the massless Dirac fermions in graphene^, there is a resurgence of interest, 
so-called supersymmetry quantum mechanics, to solve the problem in quantum mechan- 
ics in the one-dimensional Dirac equation as though in ion trapped^^ and cavity quantum 
electrodynamics^. In this scheme, the Dirac potential is a complex four-vector in which the 
bases represented in Pauli matrices. Therefore, the perturbation under this scheme is also 
represented in a vector. Specifically, we are interested to study a Dirac equation possessing 
its potential in the Lorentz force form as following (for h = 1) 



Since the study of quantum computation under the Dirac Hamiltonian is still poorly un- 
derstood, here we propose the novel technique of the generation of quantum gates by Dirac 
potential tuning in which it may be used as bandgap engineering for quantum computation. 

The change of potential and state in Dirac equation in Eq. ([1]) obeys Darboux transfor- 
mations implemented by intertwining operations on the Hamiltoniani^ii^. An important dif- 
ference from other transforms such as Fourier and Laplace transforms, Darboux transforms 
creates the new potential and state without changing their domains. Darboux transforma- 
tions has wide range applications in physics such as to find exact solutions of non-linear 
Schrodinger, sine-Gordon, and Korteweg-de Vries equations22. Darboux transformations 
for Dirac equations developed in this paper is so-called Bagrov, Baldiotti, Gitman, and 
Shamshutdinova (BEGS)- Darboux transformations^i2i in which the Darboux transform has 
an input of a unitary matrix. 

In this paper, as depicted in Fig. ([T]), we show how to generate the elementary quan- 
tum gates by the tunings of the quantum transistor potential. To do so, the transistor is 
mathematically expressed by a Dirac Hamiltonian and the tunings are done by varying the 



Dirac potential 




(1) 
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operator input, initial potential, and potential differences. In particular, one- dimensional 
Dirac equation is considered. In graphene, hexagonal stacks of carbon, one-dimensional 
Dirac equation represents the excitation of massless Dirac-fermions which are subjected to a 
uniform magnetic field perpendicular to graphene. Here, the plane wave in two-dimensional 
space is chosen for the solution of (2+l)-dimensional Dirac equation to reduce the problem 
into one dimension. Mathematical properties underlying this scheme is given in Section (jITl) . 
The tunings cause perturbation over spatial and temporal regimes of Dirac potential Bloch 
sphere. Consequently, the generation of unitary gates can be performed by controlling the 
potential differences. The exposition about this is provided in Section flllll) . 

We outline several new results of this paper: first, n-relativistic spin qubits moving in 
the external electromagnetic fields can be exploited for the generation of quantum logic 
gates. Mathematically, the quantum system is represented by a Dirac equation in which 
its potential is a Lorentz force in a complex four-vector form and the perturbation obeys 
Darboux transformation. Second, our method suggests a new paradigm of quantum gates: 
the mathematical expression for perturbing any qubit in the system and the coupling between 
two gates are well-defined, {cr^} is the coupling operator for controlling a qubit, {o"^} is 
the operator for cyclely transforming a qubit, and aao + {V{t) — /3(t))Ui is the operator for 
the target qubit. Moreover, the new type of quantum gates are proposed : U-Cyclic gate 
which can transform the qubits cyclely. 

It is conjectured that the proposal may fit with the novel dual-gate graphene field-effect 
transistor in which the bandgap is tuned by electrical field^^. 



II. MATHEMATICAL PROPERTIES 

Conventionally, the perturbation theory of a quantum system is described by the ap- 
proximation methods on Schrodinger equation^^. In this work, we propose that Darboux 
transformations on one-dimentional Dirac equations can be used to describe the theory. 

This section consists of three subsections: first. Subsection (III Ap describes the current 



problem of so-called the vector- quantum gate duality of Pauli matrices. Second, the notion 
of massless Dirac-fermion in graphene is briefly introduced in Subsection (III Bp . Third, Sub- 
section (III Cp provides the mathematical foundations of Darboux transformations on Dirac 



equations for multiple qubits. Fourth, Subsection (III Dp concerns about the mathematical 



Ay 



O O 
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FIG. 1: The proposal of voltage bias of one-qubit quantum gate generation. The figure 
shows the motivation of this paper which is to propose a novel technique of quantum 
computation transforming the qubits \initial) — )■ \ final) by tuning the Hamiltonian 
potential V + AV^ in which the is a complex four-vector representing the bandgap 
and AV^ is related to the excitation energy. The transformation is performed by Bagrov, 
Baldiotti, Gitman, and Shamshutdinova (BBGS)-Darboux transformation^^ changing 



tools for coupling of two quantum bits. 

A. Vector-quantum gates duality of Pauli matrices 

In this subsection, the duality of vector and quantum gate of Pauli matrices is discussed. 
As widely known, the Pauli matrices, {(Ji\i = 0, 1,2,3}, are the elementary quantum gates 
in the field of quantum computation^. Every quantum gates for one qubit operations can 
be decomposed into the superposition of Pauli matrices as following 



where \J(g)={g\ Identity, Hadamard, Pauli-{X, F, Z}, Phase, |} is a unitary matrix, is 
a coefficient of cTj in which it has to fulfill the normalization condition ^ |ajp = 1. 



On the other side, Pauli matrices can be used as resources for bases of a vector as 
mentioned in Ref.— . The work in this paper follows the conventions as mentioned in 
Ref.— considering the simplicity and the concept is easily understood for wider audiences. 



{V, \initial)} ^ {V + AV, \final)}. 
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The so-called complex-four vector is a vector composed by Pauli matrices as the basis vector 

3 

v = ^ Viai, (3) 

where Vi is the coefficient and cTj is the basis vector. The basis of complex four-vector, 
{(To, cTi, (T2, cTs}, obcy the exterior algebra^^ 



A a^ = ao (4a) 
(Tf, A CTj = -a., A cr«;, (4b) 

for {k, 7} = 1,2, and 3, since the basis can be divided into the scalar part, {ctq}, and the 
vector part, {cti, (72, era}. 
The vector conjugate is 

3 

(v)* = voao - ^iWi, (5) 

i=l 

therefore the Lorentz-scalar product easily reads 

3 

v.{vr = v^9,uV'' = {vor-Y,i^.)^ (6) 

i=l 

where g^^, is a Minkowski tensor. The algebraic product of two complex-four vectors is 

a. {by = a^g^^uVcTQ + ((a^6o - ^7^0) - ieai3^aj)p)a^, (7) 

where Sap-y is Levi-Civita permutation symbol. 

Eq. ([7]) has remarkable features: if the two vectors have same directions, the result is 
a scalar (or in the direction of ctq); while, if they have different directions, the result is 
a complex-four vector. Our model shows that one vector represents the physical system 
and another one represents the perturbation from environment, thus Eq. ([Tj) expresses the 
interaction between the physical system and environment. The outcomes of the interaction 
can result a scalar; it may be related to the perturbation due to the presence of Zeeman or 
Stark effect. Otherwise, it produces a vector, such as found if the perturbation is Lorentz 
force. 

In order to circumvent this discrepancy, we argue that the type of quantum gates in 
Eq. ([2]) is related to the direction of basis vector in Eq. ([2]). The correspondence can be 
emerged in the presence of the Lorentz force on the qubits. The observation of these aspects 
is provided in the following sections. 
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B. Massless Dirac-fermions in a uniform magnetic field 



Massless Dirac-fermions, massless particles which obey Dirac's (relativistic) equation, can 
be found on unbound electrons in graphene's honeycomb lattice^. The lattice is composed 
of two inequivalent sublattices labeled A and B, therefore there are two inequivalent corners 
K and K' in the graphene Brillouin zone. These two points are called Dirac points, since 
the structure of low-energy band at the points corresponds to Dirac cones^S. 

The wave function of electron which is close to K, il){r), is governed by a two dimensional 
Dirac equation 

— ivpa ■ Vip{r) = eip{f), (8) 

where vp is Fermi velocity— which has the value vp — I x 10^ m s~^ and a = {ax, cry) 
are Pauli matrices. The emergence of Dirac cones around Dirac points corresponding to 
symmetry breaking of energy band in Eq. ([8]), e, can be explained by representing the 
Hamiltonian in momentum space 

VFa-kip{k) =eip{k). (9) 

If the massless Dirac-fermions in graphene are subjected to a uniform magnetic field, there 
are degeneracy breaking of energy levels into Landau levels, which are responsible for the 
half-integer quantum Hall effect^S. Conventionally, a uniform magnetic field B which is per- 
pendicularly applied to the plane of graphene is represented by Landau gauge A=B {—y, 0). 
Hence, replacing — iV in Eq. ([S]) by — iV + ^ and for the wave function in the form 
ip{x,y) = e*'^^'0(y), we can obtain 

vp (^aydy - (k- a^^ (j){y) = e<p{y). (10) 

Eigenvalues corresponding to Landau levels in Eq. (fTOj) can be obtained by changing its 
form into the same form of atom-field interaction Hamiltonian in Jaynes-Cummings model^^ 

(Oa+ + OV-) 0(0=^10(0, (11) 

where /b = is the magnetic length, ujc = V^t^ is the cyclotron frequency of the 

Dirac-fermions, = ~5 ^ = 7^ ~ and one-dimensional harmonic-oscillator operators 

O = i= (9^ + , (12a) 



(12b) 



According to Ref.— , the A^-solutions of Eq. (fTOj) can be generated from the zero energy. 
Therefore, the eigenvectors and their eigenvalues are 



where = 0, 1,2, ... and ipN is the one-dimensional harmonic oscillator solution. The zero 
energy state of Landau levels, = 0, can only occur if quantum relativistic is taken into 
account. This particular state is useful to explain the anomalies in the quantum Hall effect. 

The above exposition shows how a particular perturbation on massless Dirac-fermions in 
graphene can explain the anomalies in the quantum Hall effect. In this paper, we introduce 
an approach in the case of the massless Dirac-fermions are subjected to any directions of an 
external electromagnetic field to perform quantum computation in graphene. 

C. Darboux transformations on quantum state and Dirac four-potential 

Darboux transformations for one-dimensional Dirac equation was suggested by Bagrov 
et al.— and Nieto et al.—. Samsonov and Shamshutdinova developed the method to show 
that it is possible to control the qubit state by the use of an external field^. Moreover, it 
was developed in Ref.~ to show the correlation between atomic inversion and Dirac potential 
in cavity quantum electrodynamics. 

Suppose that it is possible to change the potential and the state of a physical system 
obeying a Dirac equation without changing the energy of the system. Mathematically, the 
action of 'D(Uj) on a set of potential and state of a Hamiltonian is defined by V{\Ji)[N] 
{V, = {V^[A^],^[A^]} transforming the old Hamiltonian hiV)"^ = eo* -> ^(V'[l])^[l] = 
£:i\E'[l] — ... — /;,(y[A^])\E'[A^] = En'^IN]. In the context of quantum physics, it is possible 
to explain the perturbation phenomena of a quantum system due to the presence of classical 
effects using Darboux transformations: the action Darboux transformations on a quantum 
system represented by a set of potential and state, {V, \&}, cause a new set of potential and 
state of the system, {V + AV, \1/ + A\E'}. This transformation is called by Bagrov, Baldiotti, 
Gitman, and Shamshutdinova (BBGS)- Darboux transformations^ii^. Here, we modify the 



(PnAO = (0 ® It) ± (0 ® \\) , 





(13b) 
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transformation by substituting CTj by an unitary matrix, Uj, where for one qubit, Uj can 
represent a vector or Uj G {Uo= Identity gate, U{i_2,3}= Pauli-{X, Y, and Z} gates, U4= 
Hadamard gate, and U5= Phase shift gates}. 

We also modify the potential in this paper: it consists of temporal {scalar) and spatial 
{vector) terms, thus it is called Dirac four-potential. Therefore, it is mathematically defined 
by V7v(t) = Yl'j=o'^jifN{t))j- This is a complex four-vector, as defined in Eq. ([3]), in 
Minkowskian Bloch sphere and its magnitude is |VAr(t)| = r]^\fN{t))j{fN{t))i, where rj^'' is 
the metric tensor. 

Because the potential is a complex four-vector, the perturbation also follows to the exte- 
rior product or generalized cross product for arbitrary dimensions, instead of dot product. 
It means that the Dirac potential perturbation has not only magnitude, but also direction, 
which is quiet different with perturbation in Schrodinger potential consisting the magnitude 
only. The consequences of this choice are shown in the next section. 

^^j^ 2x2 4x4 8x8 16x16 

One can consider to use the generalized Pauli-Dirac matrices, Oj = (^,^^,^n~^, ,...}, 
where the overbraces denote the dimensions of the matrix, as given in Ref.— , if the space- 
time dimensions are higher than 4 and also to expand the Dirac equation operating on 
n-qubits. 

The nature of a Dirac equation admits the simulation for two qubits. As shown in Ref.—, 
by assuming that the particle is massive^, m ^ 0, and in the rest frame of particle, A;^ = 
(m, 0), the Dirac equation implies that the set of a quantum bit {|00), |01)} is belonged to the 
positive energy, ip^^\ ^iiid the set of quantum state {|10), |11)} corresponds to the negative 
energy, ip^~\ if the solutions of the Dirac equation are plane wave, 

this work, we generalize the notion into n-qubits and also the potential is modified into higher 
dimensions. Therefore, a modified Dirac equation representing from one qubit, |({0, 1}^)), 
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to n-qubits, |({0, l}")), is 

V 



3 



V 

3 



(=0 



+ =^I({0'1P)) (14) 



y 

3 
1=0 

; i 

Hp'^a, + J]/,p,)l({o,ir)) =^l({o,ir)) 

1=0 

i i 

V 

3 
1=0 

Eq. (fT4|) has several possibilities of interpretations. According to the original Dirac 
theory, n-qubits would need ^ Dirac-fermions. However, in conformity with the model of 
massless Dirac-fermions in graphene as shown in Subsection fill Bl) . n-qubits correspond to 
n-massless Dirac-fermions. 

Due to the perturbation, it changes the potential Vo(t) = Xlj ^i(/o(^))i ^ ^i(^) = 
~^ ■■■ ~^ VN{t) = J2j ^jifN(t))j obeying the intertwining operation 

£{V,)h{Vo) = h{V^)C{lJ,) ^ C{lJ,)h{V,) = h{V2)C{lJ,) ^ ... ^ C{lJ,)h{VN-i) = h{V^)C{lJ,). 

(15) 

The chain of relations in Eq. ( IT5|) means that the new potential, from Vi to Vat, can be 
generated under intertwining operations from the initial potential (Vq). The intertwining 
operator is 

£(U.) = | + i5(U.), (16) 
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where 



i3(U,) = a,{t)ao + {Vr,^i{t) - A(t))U,. (17) 

In Eq. (fT7|) . Uj is a vector represented in matrix form. These also lead the linear transfor- 
mation of state2i, for one qubit, ^ ^ ^[1] = £(Ui)^ = £(Ui)^[l] ... *[iV] = 
£(UO^[iV-l]. 

It is clear that the terms {oi, represent the perturbation terms. To conclude about this 
method, in other words, the change of Dirac potential under BBGS-Darboux transformation 
obeys intertwining operations, while the change of state follows linear transformations^. 



D. Control and cyclic operators 

In this subsection, we introduce the two operators which will be extensively used in 
constructing quantum circuit. These operators satisfy idempotent operation, = A. First, 



control operators, {a 



(±) 

(0,3) 



(±) 

(co) 



}. Second, cyclic operators, {a 



(±) 

(1,2) 



(±) 
(cy) 



ai:Lio'2 



which are widely known as matrix representations of the raising and lowering operators^^. 
Let us define: 



co{±} 



crltl (S)A+ all ® B, 



'{CO) 



'{co) 



:i8a) 



cy{±} 



a 



(cy)' 



;i8b) 



It is shown below that these formalisms are useful to construct any unitary matrices and 
to reduce a lot of space to write multiple qubit operations. From these operators, one can 
obtain the following gates for two qubits: 



1. SWAP-gate = aj+j ® a'itl + a^"^ 



'{co) 



'{CO) 



{co) 

(+) 

(cy) ^ , (cy) 
cy{±} 



^(5 © ^(U + © ^ 

' co{±} 

/l o\ 
10 
10 
\0 1/ 



cr 
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2. Full SWAP-gate = 



a 



(-) 

(CO) 



(+) 

(CO) 



(cy) 



^(cy) + ^(cj;) ® ^{cy) 



{cy) ^ (co) (cy) 



a 



(+) 

(co) 



0" 



(+) 



(cy) ©/^(U 

cy{±} 



^ 

cy{±} 

/o l\ 

10 
10 

\i oy 

It is named by Full SWAP-gate, because it transforms |00) o |11) and |01) ^ |10). 

Let us consider the use of control- and cyc/zc-operators for constructing the quantum 
mtes. 



Controlled-U gate = cr|^^^ ® ao + cr| (g) Ui 

co{±} 



(19) 














1 




























U-Cyclic gate = ctq ® crj^y) + crg^ 

= ^0 U,. 

cy{±,T} 



(20) 



If Uj is (Ti and by choosing Eq. ([20]) in the form do ® a^(^J) + Uj ® a^^J) = ctq ®cy{T} 
the gate is called clockwise-cyclic-ga.te (CC-gate), because it has the form 














1 














1 














1 






and it transforms 



|00) ^ |01) 
|01) -> |10) 
|10)->|11) 
|11) ^ |00) . 



(22) 



The counterclockwise-cyclic-gate (CCC-gate) can be performed by the substitution of 



Uj by (Ti and Eq. (120|) in the form ao cr^ 
as following 

/ 



{+) 



{-) 



and it changes 






1 














1 














1 


1 












Co © Uj. It has the form 

cy{±} 



(23) 



|00) ^ 111) 
|11)->|10) 
|10)-^|01) 
|01) ^ |00) . 



(24) 



In other words, (CC) = (CCC) = 1 is to perform one round rotation as shown in Fig. 
([2]). Interestingly, one can also find the following relations 



(CC) = (CCC)T = (CCC)"^ 
(CCC) = (CC)^ = (CC)"\ 



(25a) 
(25b) 



which are the properties of unitary matrices, besides the facts that |det(CC) | = |det(CCC) |=1. 

The definitions in Eq. (fTS!) can also be used to represent the standard sets of universal 
gates of one qubits 



Hadamard gate = (a';tl - cril + cr'^tl + 



^2 V (™) (^2^) ^"y'^ J ' 

(co) ^ (co) ' 



Phase gate = a^^^l + ia^ 



(26a) 
(26b) 



13 











111) |01) 




111) 


101) 


^ 110) ^ 






10) <^ 


2a Illustration for clockwise- cyclic 


2b Illustration for 



operations of two qubits. 



counterclockwise- cyclic operations of 
two qubits. 



FIG. 2: The illustration of how the use of U-Cyclic for quantum computation is. If Uj is 
(Ti, one can obtain clockwise- cyclic and counterclockwise- cyclic operations of two qubits. 



TT 

4 



(26c) 



Controlled 



NOT 



and for two qubits 

Controlled-NOT gate = 0^^^- (27) 

co(±) 

The application of the approach in this Subsection to the rest of discussion is given in 
Subsection fUTTDl) . 



III. QUANTUM GATES GENERATION BY A QUANTUM TRANSISTOR 

The relativistic electron theory of Dirac successfully explains the quantum phenomena 
such as the g factor and the proper fine structure of Zeeman effect correctly^. Several 
efforts have been achieved for a better understanding of Dirac theory by modifying the 
equation. The famous one was presented by Foldy and Wouthuysen which was called by 
Foldy-Wouthuysen transformation afterwards^^; they showed a kind of canonical transforma- 
tion in which can connect between Dirac (or relativistic quantum) and Pauli (non-relativistic 
quantum) theory and to generate Zitterhewegung. Nevertheless, it is still elusive what kind of 
transformation on Dirac equation can be used to implement quantum computation on Dirac 
theory. Motivated by supersymmetric quantum mechanics, we shall present the technique 
to exploit the one-dimensional Dirac equation into the field of quantum computation. 
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Suppose that it is necessary to transform a Dirac equation possessing a classical pertur- 
bation as following 

hoiV) hiiV+AV) 

\-ze4, + v)\{o, in. = ^o|{o, in, ^\-ze''d, + v + Av)\{o, in f = £i|{o, in / 

(28) 

where {V, |{0, is a initial set of potential and qubit and {V^ + AV^, |{0, l}")f} is a final 

set of potential and qubit. The matrix 6'^ is a generalized form of Pauli-Dirac matrix as 
explained in Ref.— , 6*^ = {cr^, 7^^, p'^, ...}. Therefore, Eq. fl28l) provides a correlation between 
Dirac potential and n-qubits quantum computation if one-fold Darboux transformation on 
the equation is necessary to be defined. 

Below, we present that this kind of transformation can be realized by the action of BBGS- 
Darboux transformation as introduced in Subsection fill CD on the equation possessing the 
complex four-vector potential. It is shown below that quantum gates can be generated 
due to the perturbation of the potential of n-relativistic spin qubit moving on classical 
electromagnetic field. The exposition for a single qubit, two qubits, and followed by a 
generalization to multi-qubit systems are given; also, another way to construct a quantum 
circuit, i.e., in the case of coupling one relativistic qubit, is presented. 

A. One qubit quantum gates generation by a quantum transistor 

In this scheme, initially the particle is at rest and has the Dirac potential Vq and in 
quantum state |{0, This set of Dirac potential and state, {Vq, |{0, belongs to 

the following one-dimensional stationary Dirac equation 

ho\{0,lY), = eo\{0,iy)i, (29) 

where ho = (ia^j^+Voit)), Vo{t) = Zlj crj(/o(t))j, and Eq is a constant. 

The unitary transformation on the initial quantum state resulting the new quantum state 
is achieved if the relativistic qubit moves in the external electromagnetic fields so that the 
system has a set of new potential and state, {Vq + AV, |{0, 1}^) f}, which belongs to the new 
Dirac Hamiltonian 

h\{0,iy)f = e,\{OAy)f, (30) 
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where = {iaz-^+Vf{t)), Vi{t)=Vo + Al^=^^. and ei is a constant. 

Let us assume that the potential perturbation, AV, is related to a relativistic qubit with 
charge in external electromagnetic fields, hence the new potential is in the following form 



Vi = Vo + AV 
P_ 



Vo + (31a) 



3 

where p=J2 (p)"'^a reads^ 

a=0 



= qv^{d''A^ -df^A") (32a) 

and where q= charge of the qubit, = {vq, —Vi, — f2, — "^3) is the four-velocity of the qubit, 
= {A^ , A^ , A"^ , A^) is the four-potential of the external electromagnetic fields, and F"^ 
is the antisymmetric field-strength tensor. In the complex-four vector representation, the 
Lorentz force is in the following forn>2^i2^ 



p= -q 



v.OA + DA.v 



2 

= Re{v.DA) (33a) 

Therefore, the representation of Eq. ([T]) in a complex four-vector form is 

3 

^^^0^9, + J2 fkOk)m 1}")^ = ^o|{0, 1}").. (34) 

fc=0 

For the case of a single qubit, the perturbation is represented by the exterior product 
between two complex four-vector represented by 2 x 2 matrix as following, 

AV = -la, A Ui, (35) 

where the eft. belongs to the direction of relativistic spin qubit and Uj is related to the 
direction of external electromagnetic fields. Therefore, the direction of external classical 
electromagnetic fields is an input of the perturbing potential. 

The proposed scheme covers potential initialization and switching. In the context of 
dual-gate graphene field-effect transistor, the initial Dirac potential may be related to the 
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perpendicular electric field generating the bandgap^ and the switching Dirac potential cor- 
responds to the energy transition of the electrons needed for the excitation from the valence 
band to the conduction band. It is shown in the theorem below that the process of potential 
initialization and switching is unique for every quantum gate generation. Following exposes 
the detailed explanation of the scheme. 

Theorem 1. Suppose the initial potential is Vo{t) = o"i(/o(^))j- The one fold BBGS- 
Darhoux transformations on the equation / [gPj) at which the final potential is Vi{t) = 
Vo{t) + AV,where AV = — icrt. A Ui, is a suffice condition for >C({Uj}) = {Uj}. (Vo(t))o 
and (Vi(t))o are the vector variables. 
Proof. Consider 

£(Ui)/io^ = hjC{V^)^. (36) 

One can find 

MVoit) - Vo{t)) + t[lJ„ a,]{Vo{t) - m - 1) = (37a) 

la^jSiit) - V,{t) - o^i) + aMV,{t) - V,{t)) + {V^{t) - P,{t)){V,{t) - (Vo(t)U, 

+^c^.((/3o(t) - Vo{t) - ai{t)) + aoK)(t). 

(37b) 

Let us examine the constraints of how to generate Uj gate. By setting Ui up to 1, the 
resume for Uj can be obtained as shown in Tabs. ([T]), ([Tll), (IIIip . and (lIVp . respectively. 

The illustrations show the initial potential preparations [blue dashed line) followed by 
the tuning of final potential {red dashed line) to generate a certain quantum gates. The 
tunings can be accomplished by several manners: first., temporal displacement of the initial 
potential, second., displacement on the potential magnitude, and, third, by changing the form 
of the initial potential function. Under the scheme, some light is cast on the question of how 
a quantum transistor achieves quantum computation by an examination of the correlation 
between Dirac potentials and quantum gate generation. 

As shown in Tab. ([T]), the IDENTITY-gate can be realized by preparing the initial po- 
tential in the magnitude of parabolic form. The structure of potential function is transformed 
from time-dependency to space-time dependencies due to the magnitude displacement of the 
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A 




FIG. 3: The illustration of the possibility for the implementation in novel dual-gate 
graphene FETs^^: a, A non-zero bandgap A and a the Fermi energy shift, Ep-, are induced 
by the electric fields, b, The perturbation on Dirac potential, AV^, excites the electron 
from the valence band to the conduction band and performs the quantum computation. 

potential. It is similar to the case of quantum state transition due to carrier-photon scat- 
tering event in intraband of semiconductor-^: the relativistic spin qubit absorbs the light 
during the moving. By assuming the photon momentum to be zero, the transition on the 
band diagrams are always vertical. As illustrated in Fig. ([3]), in the case of dual-gate 
graphene field-effect transistor, the parabolic Dirac potential yields the bandgap, then the 
electron is excited to the conduction band due to the tuning of electromagnetic field. 

The NOT-gate in 2;-basis can be performed by temporal displacement of the initial to 
the final potential in the parabolic form. The resume of this scheme is given in Tab. ( IIT|) . 

Tab. f lllll) provides the scheme to generate CTy-Pauli matrix. It can be realized by sharp- 
ening the initial parabolic potential followed by temporal and magnitude displacements of 
the potential. 

Furthermore, the (Jz-Pauli matrix which is NOT-gate of phase state can be accomplished 
by setting the initial potential up in constant function, followed by tuning the final potential 
in parabolic form. 

The other type of quantum gates can be determined by similar manners. 
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Quantum gate 


Input 


Initial potential 


Final potential 


Ui 
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= t((T2 - ifJs) 
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TABLE II: Resume of Ui 



B. Multiple qubit quantum gates generation by a quantum transistor 

Remarkable fact about 7-matrices is they have similar mathematical structures with Pauli 
matrices^. Quantum gates for two qubits can also be decomposed into 7-matrices2^. 

To demonstrate how the proposed principle works in the case of two qubits, one can use 
the 7-matrix— to describe the complex-four vector as following 

3 

/i=0 

where is a complex-four vector represented in 7-matrix and symbolized by Feynman slash 
notation. 

It is also clear that the product of two vectors, a^fe'^, in the representation of 7-matrix 
obeys complex four-vector since the 7-matrices satisfy Clifford algebra 

{7^ 7"} = 7''7" + 7^7'' = 27/^^^ (39) 
The initial one-dimensional Dirac equation of two relativistic spin qubits at rest is rep- 
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Quantum gate 


Input 


Initial potential 


Final potential 


U2 


0-2 


Vo{t) 
= -t{ai + ia-i) 


= -{t + l)(Ti - itas 



TABLE III: Resume of U2 . 



resented by 

(^3^ + Vo^ |{0, = eo|{0, 1}^)^, (40) 

where {Vq, |{0, is a set of initial potential and quantum state of the system consisting 

two qubits. 

The perturbation on the system occurs if the two relativistic spin qubits moving in clas- 
sical electromagnetic field, therefore the quantum system obeys the new one-dimensional 
Dirac equation 

(^73^ + Vo + AV^ |{0, = £i|{0, 1}^)^, (41) 

where {Vq + AV, |{0, 1}^)/} is a set of final potential and quantum state of the system. 
Similar to in the case of one qubit in Eq. (13T!) . here the perturbation term corresponds to 
Lorentz force, therefore q represents the effective charge of the two relativistic spin qubits. 

The change of potential and state in the case of two qubits follows the Theorem 1, 
nevertheless the intertwining operator in Eq. (fTTI) is substituted by 

i3({C(U,), SWAP}) = a,(t)70 + (Viv-i(t) - A(t)){C(UO, SWAP}, (42) 
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TABLE IV: Resume of U3 . 



where C(Uj) is Controlled-Uj gate, and the perturbation on Dirac Hamiltonian is defined 
by 

AV = -27l A {C(U,), SWAP}. (43) 

Similar to the case of a single qubit, the term {C(Ui), SWAP} in Eq. (l43l) corresponds 
to the direction of an external perturbation, while 7^ belongs to the direction of the two 
relativistic spin qubits. This method converts the "vector basis" aspect as found in Eq. ( I43p 
into "quantum gates" aspect of 7-matrices, since the initial Dirac potential is vanished by 
the perturbation term due to the transformation. The scheme of this method is illustrated 
in Fig. (H. 

By similar manner, one can obtain for n-qubits, since the vector basis and quantum gates 
are easily acquired by Pauli-Dirac matrix generator—. The multiple qubits quantum com- 
putation on Dirac equation by utilizing Lorentz force can be realized by simply substitution 
of Pauli matrices in the Theorem 1 by 6'^-matrices. 

To generalize our proposal, the following set of equations show the multiple qubits quan- 
tum computation on Dirac equations utilizing Lorentz force which is transforming the initial 
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Ay = -i7l{C(U,),SWAP} 



Cho — hiC I 



Vq + ^V 



£({C(U,),SWAP})* = 

vector basis 



->{C(U,,),SWAP}|{0,l}2), = |{0,lf 



quantum gate 



FIG. 4: The scheme of methods for two qubits. The calculation inside the red box 
represents the main method: Darboux transformations on one-dimensional Dirac equation 
causes the transformation a set of potential and state from the old into the new one, 
{Vq, |{0, lY)i} {Vq + AV, |{0, 1}^)/}- However, if the potential perturbation term is 
fixed in a certain formulation related to Lorentz force form, as written inside the yellow 
box, the transformation operator of the quantum state is changed into the form of 

quantum gates. 



set of Dirac equations as given in Eq. (fT4|) into the final set of new Dirac equations 



{-i$+fhm, m)^ = eom, ^ {-i$+fmo, i}^)>/ = eii({o, iv))s 

{-i$+fb)m,m)i=eo\miY)h ^ M^+z/oiacip))^ =£ii({o,i}2))^. 

i AV = -i-fi A {C(U) ., SWAP} I 

M^ + //&)I({0,1}^)>. =eo|({0,iy')>, -> (-i^ + //i)|({0,l}=*))/ =£i|({0,l}3))^. 

4, AV = -i7r3A{TOFFOLI,CC(U)i,...} I 

+ /o)l({0, l}")>i = eol({0, ^ + //i)[{{0, 1}"))^ = £i|({0, 1}"))^ 

4. Ay = -163 A {Controlled"-i-(U) ., ...} i 
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C. Coupling of quantum gates 



According to tlie previous sections, for instances, the ControUed^-NOT can be obtained 
by tlie following expressions 



Controlled^-NOT = CNOT = Uq Ui (45a) 

co{±} 

Controlled^-NOT = TOFFOLI = Uq Uq Ui (45b) 

co{±} co{±} 
N-1 

, V 

Controlled^-NOT = Uo - Uo Ui (45c) 

co{±} co{±} 

Then, for NOT-Cyclic^, one can find that 



NOT-Cyclic^ = Ui Uq (46a) 

cj/{±} 

NOT-Cyclic^ = Ui Uo Uo (46b) 

CJ/{±} cj/{±} 

. ^ s 

NOT-Cyclic^ = Ui Uo ...Uo (46c) 

cy{±} cj/{±} 

D. Physical interpretation of control-, cyclic-, and target-operators 

In this section, we present a new interpretation of the control-, cyclic-, and target- 
operators in quantum gates based upon the explanation in the previous sections. The 
illustrations are given in Fig. ([5]) for the quantum gates of two qubits. 

Under this scheme, target operator acts on each qubit and is coupled by {control, cyclic}- 
operators for constructing n-qubits quantum gates. 

It is shown in Fig. flSa|) that when the control operator is turned on, the circuit is in 
Controlled-NOT form and it is mathematically expressed by Eq. (IT^ . 

As can be seen in Fig. (I5bp that the coupling operator is changeable, then, the turning 
cyclic operator on causes the gate in NOT-Cyclic form where it is represented by Eq. (|20|) . 
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5a Illustration for Controlled-NOT. 




5b Illustration for NOT-Cyclic. 

FIG. 5: The illustration of our mathematical proposal for describing the process of 
quantum computation due to the perturbation of the target qubit by an external force. 
The target qubit performs Unitary gate due to the perturbation. One can realize 
Controlled-NOT by setting the control operator, while for NOT-Cyclic by cyclic 

operator on the coupling. 

Another interesting feature of this scheme is that, for a single qubit, both positive and 
negative energy, {tp^^\ip^~^}, can be used one-qubit quantum computation. As has been 
discussed in Sec. (Jill), the degeneracy is broken for the case of multiple qubits: the positive 
(negative) energy, ip^^^ {^^~^), is the signature of |0) (|1)) state of the first qubit. Therefore, 
swapping the quantum state of multiple qubits is simultaneously interchange the energy of 
the system under the action of {SWAP, Full SWAP} gate which is realized by the setting 
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of cyclic-opeiatoT on the coupling, while the chose of control-operator ensures the energy 
conservation. 

IV. CONCLUSION 

In this paper, we propose a modification form of a Dirac equation by allowing the potential 
in complex four-vector form representing the classical perturbation on a quantum system. 
The notion of this work is similar as proposed in Ref.— : the modified Schrodinger equation 
form can be used to represent nonlinear phenomena of a classical system such as water wave, 
by allowing the potential and also its perturbation are in the classical form. 

We have shown that ra-relativistic spin qubit in external classical electromagnetic fields 
can be used for resource of quantum transistor. The mathematical explanation behind the 
proposal is on the action of Darboux transformation on one-dimensional Dirac equation at 
which the potential is a complex Dirac four-potential utilizing the Lorentz force as perturbing 
potential. These mathematical methods take advantage of the vector- quantum gates duality 
of Pauli matrices. 

It is also shown that our proposal can closely describe the resources of quantum tran- 
sistor required for implementing quantum computation: every n-qubits require n-quantum 
transistors. In this scheme, the transistor is prepared with certain initial potential function, 
Vinitiai, then the final potential function, Vfinai, is obtained to generate certain quantum 
gate. 

However, the proposed scheme as presented is aimed as a primitive approach of a quantum 
transistor based upon Dirac equations. The next challenge is how to implement it into a more 
specific and complex physical system. The correlation between a specific physical quantity, 
as though quantum Rabi oscillations in cavity quantum electrodynamics and bandgap in 
various types of graphene transistor, and the change of qubit state due to the tuning of the 
direction of the external electromagnetic field, shall be explored. 

The proposal may be also very useful for explaining the quantum system used not only 
the magnitude but also direction of a physical quantity to encode the qubits, as found in 
superconducting flux qubits, where the two distinguishable configurations are obtained by 
induced clockwise or anti clockwise current^. 

The proposed formalism opens the possibility to involve the relativity into the further 
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study of quantum computation and information, coincidentally with an emerging assertion 
that closed timehke curve can solve classical problems in computation such as NP-complete 
problems^"—. 

For the next work, it is also very interesting to consider if quantum circuit is constructed 
under the Dirac equations for ri-qubits instead of conventional way in Schrodinger represen- 
tations. The similar evaluation as provided in Ref.— may give noteworthy results since this 
kind of construction may influence the depth and space resource of a quantum circuit. In 
this work, we assume that intertwining operation runs only the potential of relativistic spin 
qubits. The computation of a qubit still runs under linear transformation. If the intertwin- 
ing operation runs both on the transformation of potential and qubit, this may contribute 
into the theory of the topological quantum computing^ on a massless Dirac-fermion. 
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